Lehmer's interesting series 



Freeman J. Dyson 

Institute for Advanced Study, Princeton, New Jersey, U.S.A^ 

Norman E. Frankel 

Department of Physics, University of Melbourne, Victoria, Australi^ 

M. Lawrence Glasser 

Department of Physics, Clarkson University, Potsdam, New York, t/&4@ 

(Dated: March 29, 2011) 

Abstract 

The series 



Sk(z) = ]T 



— m k z m 



j (2m\ 
=1 \m) 

is evaluated in a non-recursive and closed process and it can be analytically continued beyond 
its domain of convergence < \z\ < 4 for k = 0, 1, 2, • • • . From this we provide a firm basis for 
Lehmer's observation that tt emerges from the limiting behavior of S&(2) as k — > oo. 
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I. INTRODUCTION 



In the paper jlj marking the 60th anniversary of his first contribution to the American 
Mathematical Monthly, D. H. Lehmer studied two classes of interesting infinite series; in- 
teresting in the sense that the terms and sum are known explicitly. Lehmer evaluated the 
series of class II, whose origin he attributes to L. Comtet 0, namely 



S k (z) = J2 



m k z m 



m=l 



/2m\ 



(1) 



for non-negative integer k. He did so by manipulating the Taylor expansion for the arcsine 
function, so arriving at a general formula for Eq. (JTJ in terms of recursively defined polynomi- 
als. As we show below, the domain of convergence for Eq. (|Tj_ is \z\ < 4. Among publications 
stimulated by Ref. [l| are those of Borwein and Girgensohn 3| and of Mathar [4J. Also there 
are recent papers [5Hl4j that contain other relevant references. 

Lehmer's class II series can also be expressed as generalised hypergeometric functions 
and this has been exploited by Krupnikov and Kolbig [15J to construct a table of values of 
p+iF p for unit argument and rational parameters. Lehmer's result, 



S 2 (2) = 11 



|tt, 



(2) 



illustrates the fascinating feature that 5^(2) = R\{k) + R2(k)iT, where Rj{k) is rational and, 
empirically, the ratio Ri(k)/R 2 (k) is an approximation to tt which improves as k increases. 
In fact, as we shall illustrate below, in this way Sfc(2) appears to give tt to roughly k decimal 
places. Our aims in this paper are to express Sk(z) for \z\ < 4 in non-recursive form and to 
investigate the indicated approximation to tt. 



II. CALCULATION 

We begin with the observation that — C"^ 1 = B(m,m + 1), where B is Euler's beta 
function and denotes the binomial coefficient. Hence 

Sk(z)= / ^m w [ri(l-t)] B . (3) 

Jo 1 m=l 



Next, Euler's famous 1775 discovery [161 ]. embodying the Eulerian triangular numbers, 

oo p n 

™ P X m = (-l) m+n m p X n (1 - X)- n ~\ (4) 

m=l n=l m=l 

gives 

It can be shown that 

fUt [zt(i-t)r _ r(n) , 



2 



where X = zj (4 — z) so 



fc+1 n 



m+n 



n=l m=l 

This gives, for example, Si (2) = n + 3. 



£W vnk+1 p / 1,1. y\ 

r( n + i) 2-Pll-5> n > n + 2' 



(7) 



Now using Eq. [7.3.6(1)] in Ref. [Uj, i.e. 



2 Fi(-i,n;n 



|; -1) = 2- v^F r(n + 1) { [r(n±i)]- 2 + | [r(f + I)]" 2 } , (8) 



n=l 



where 



E(k,n) 



[r{*±)Y 2[r(^)] 

(_l)n ™ 



2 [ ' 



n! 



£(-irc, 



n m fc+1 



(9) 



(10) 



m=l 



are Stirling numbers of the second kind. The right hand side of Eq. is manifestly of the 
form R\{k) + R 2 (k)n and leads to Lehmer's insightful observation, which we express more 
transparently in the form 



where 



1- '41 = 2. 



ti{k) = t u {k)+t 12 {k) and t 2 {k)=t 2l {k)+t 22 {k) 

E(k,2n) (2n)\ (2n - 1)! 



t n (k) = 



n=l 



[r (n + m 



^ E(k,2n + 1) [(2n + l)\} z 



n=0 



n=l 



(ii) 

(12) 
(13) 

(14) 

(15) 

(16) 



n=0 



For large summand n it is clear that the sums t^- are dominated by their "tails", but for 
large n 



(2n-l)! ^ l_ c2n 
2 " 



[r(n + | 



1 [ (2n + 1)! 

2 |_2r(n + §)_ 



->■ Cf n . 



Hence, as — > oo, 2t n (/c) — > t 2 i(k) and 2t 12 {k) — >■ t 22 (k). Thereby Lehmer's assertion is 
demonstrated. 
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For z = —2, we have 
fe+i 



S,<-2) = £<-!)■ n< ^ 



?1=1 



I ^ (±) p 2* 



i r(0 



Next, generalising Eq. (jSJ) (Glasser unpublished), 
2 -Pi(-|, n;n + |; -2) 

"-' ' -l) k T{k + \ 



1 1 



2 / /-w 



1 /z + l 



fc+1 



x < \/ z sin 



z + 1 



if ('-l)! 
2 ^ \ : + 1 



1=1 



2/ n nl 



1 + n(z 



+ 1) Sin "(\/7TT) ^i(i-n.i;2;*?) 

ylA(i)»rffl (1-n), ^ + l^- m 

fc=0 



tr ( 2 )*(§)^ 

The hypergeometric function above is a Jacobi polynomial [l 

iV! „(!,-*-§) 



2 F 1 (-iV,i;2^ 



(2) 



A? 



A? 



(1-2*). 



Therefore, we have our principal closed form result 
fc+i 



n=l ^ 







X 



n— 1 / 1 



E(k,n) 



p+1 



(17) 



(18) 



(19) 



4-, \2j 2^(|)^4 



From Eq. fl20|) . /or rational z 



X1 m k z m 



f2m 
m=l m 



R 3 (k) + R A (k)J- sin" 1 -V 



2"' 



(20) 



(21) 



where, again, Rj(k) is a rational number. 



III. DISCUSSION 



As pointed out in the introduction, the radius of convergence for S] C (z) is 4. Furthermore, 
since for large m the asymptotic behavior of the m— th term is m k+l l 2 (z/A) n \ the series 



4 



converges nowhere on its circle of convergence for positive k. However Eq. ( 120]) shows that 
Sk(z) is analytic on its two-sheeted Riemann surface joined along a branch cut emanating 
from z = 4. The corresponding values on the two sheets differ only by the sign chosen for 

2 we have 5^(2) 



\/4 — z in Eq ( 1201) . For z = 2 we have <Sfc(2) = Ri(k) + R,2{k)7r where the rational number 
Rj(k) has denominator 1 or 2. We show values of the ratio, Ri(k)/R2(k) and of S k (2). In 
Table [II] (and continued in Table 11X11 we list the ratio for k = 1, 65 to 65 decimal places. 
The last line of Table [TTT1 contains the first 65 decimal digits of the value of tt. For k = 100, 
96 places of tt are reproduced. In Table HVl we display the first thirty complete expressions 
for S&(2). As has been demonstrated in the previous section, the limit k — » oo associated 
with S k (2) is indeed tt. The rate of convergence is derived in the next section. 

The only previous non-recursive evaluation of an S k series appears to be that of Borwein 
and Girgensohn [3] for their b 2 (k) = Sjt(l) found by a different procedure to ours. By 
comparing their formula with Eq. |2"U|) we obtain the intriguing identity 



fc+i 

E 

n=l 



-1) 



n+k+l 



7V. 



3" 



E{k,n)C 2 n 



2n 



n-l 



fc+1 



n=l 



n 



--E 

P=0 



x ^ 3 P 

^(2p+l)C* p 

(i) 4P+ 1 
1 (p + 1)! 



2tt 



p=0 



n-l 



Eh 



r(() 



7T 



)*4< 



27 



(22) 



Note, that by equating Eqs. (EJ) and (120]) for 2 = 2, one obtains a similar identity. 

As we wrote at the outset, we were stimulated by Lehmer's paper jlj, and especially 
his proposal that S k {2) is intimately related to tt. It is true that this way of computing tt 
by taking k — > oo is not as productive as others, such as Ramanujan's modular equation 
approach as given in a nice review [HI]. Nonetheless we believe that the approach is a 
fascinating contribution to the lure and lore of tt; one which we have found both motivating 
and rewarding. In his honour we have named Um.k^ooRi{k) / R2{k) = tt, Lehmer's limit. 



IV. THE ERROR 

To derive the error to the k th approximation for tt, E{k), we introduce the difference 
function, Dk(2). Eq. ( |20l) showed that the function Sk(z) is analytic in the entire complex 
plane with a cut along the real axis from 4 to oo. Furthermore, the function can be continued 
analytically across the cut onto other sheets where it is still analytic. It is easy to show that 
the sum, Sfc(2) = R 2 {k) tt + Ri(k), is the value of the function at the point z = 2 in the first 
sheet, while the difference, Dk(2) = R 2 (k) tt — Ri(k), is the value of the same function at 
the point z = 2 on the second sheet. 

These two functions define the error, 

1 ) R 2 (k) S k (2) + D k (2) ~ S k {2) • 1 ) 

We now introduce their generating functions, whose definitions are 

oo k oo k 

G(f) = 5>*(2) n ; H{t) = Y,D k {2)-, (24) 

fc=0 " ° fc=0 ' * 
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where 



and 



G(t) 



H(t) 



[sii 



s sin 



(1-S 2 )5 

s [cos _1 (s)] 



+ 



s 2 y 



(25) 



(26) 



with s = ~^ e2 - 

It is remarkable that the two functions look similar but have very different behavior. G(t) 
has its closest singularity as s — 1, with t = ln(2), while the closest singularity of H{t) is 
at s = —1, with t = ln(2) + 2iri and ln(2) — 2ni. For large /c-values, only these leading 
singularities need be considered. 

Adding Eqs. ( 1251) and ( |26i) gives 



i 

2[2" 



, 3 ' 

,12 



and therefore 



A;! 

47TZ 



C2 



(2 



(27) 



(28) 



After shifting the singularity with t = ln(2) + x, and introducing the function g(x) which is 
finite at the origin, 



x 



with a few of the coefficients given by 



C2 



3=0 



X' 



(29) 



3 


c j 





1 


1 


-1/4 


2 


-1 / 32 


3 


5 / 384 


4 


7 / 10240 


5 


-19 / 40960 


6 


-869 / 61931520 


7 


715 / 49545216 


8 


2339 / 7927234560 


9 


-200821 / 475634073600 


10 


-12863 / 1993133260800 



we have 



'-x)» {ln(2) +x} 



k+l 



-dx . 



(30) 
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The contour can be moved to surround the positive x-axis (Hankel contour) and, accom- 
modating the jump across the cut, Eq. (!30|) becomes 



R2 



3=0 



where 



\k+l 



dx. 



(31) 



(32) 



'o (1 + x) 

The integral, Ij, can be evaluated with Euler's beta-function and gamma-function reflection 
formula, providing the result 



#) = z E c j {in(2)r fc - 



3=0 



r(fc + f-j) 
r(l-i) 



(33) 



where C,- 



Similarly, D 2 {k) can be obtained, by first observing that Eq. (|2"7|) for could have 

been obtained also by taking only the first term for G(t) in Eq. (1251) . with sin _1 (l) = |. 
Thus we have directly from Eq. (1261) 



[2-e* 



2 ' 
_> 



Everything goes through in exactly the same way as for R 2 (k), providing the result 



D 2 (k) 



-2tt3R 



.i=0 V2 



(34) 



(35) 



Our final result for the error, E(k), in —expansion form can be obtained readily from 
Eqs. (E3J and (E5D, as 



E(k) 



> 7T 



3? 



E^ F i {lB(2) + 27ri}- 



(2fc+3-2n) 



where Fj = — (2j — 3)!! Cj, (—3)!! = —1, and (—1)!! = 1. A few of the coefficients are 



(35a) 



J 


F 3 





1 


1 


1/4 


2 


1 / 32 


3 


-5 / 128 


4 


-21 / 2048 


5 


399 / 8192 


6 


869 / 65536 


7 


-39325 / 262144 


8 


-334477 / 8388608 


9 


28717403 / 33554432 


10 


59697183 / 268435456 
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Note, no Stirling's asymptotic expansion is required in obtaining Eq. (135aj) . 

The order of magnitude of E(k) can be seen by looking at just the j = term of Eq. (I35al) . 
which is 



where 



tan 



and 



E (k) 



8tt 



ln(2)7 



Q 



k+i 



COS 



k + 



4609 ■•• rad = 83.7047- ■■ deg, 



(356) 



Q 



9.1197' 



2vr 



8tt 



As ~ 1, the order of magnitude of the error is 



E(k) = O 



2tt \ 



-fc' 



(e.g. 10~ 96 for k = 100) 



(35c) 



To assess the overall accuracy of the expansion of E(k) given by Eq. (135a[) . we have 
compared values found from it with just the first seven terms ( j = — 6) with the exact 
numerical result for it — -f^y given in the previous section. That comparison is shown in 
Table [fl For large values of k, the agreement is excellent and that remains very good for 
lower A;- values. The question of the convergence of the expansion for E(k) in Eq. ( 135a|) . and 
the continuing considerations this question engenders, are discussed in the next section. 



cj oc [2tt] 



V. CONSIDERATIONS 

The coefficients Cj behave in an interesting way. The even and odd coefficients behave 
differently. They vary as 

for even j 
y/j for odd j 

After the first few, each odd-j coefficient is larger than the preceding even-j one. This occurs 
because the generating function for the series has singularities at plus and minus 2m which 
alternately give terms that add and subtract. Similar things occur with the Lehmer errors, 
E(k), but the periodicity is not so obvious as the phase shift is not exactly |. 

It is clear now that that the series for B,2{k) is a convergent one, while that for D^(2) 
is divergent, albeit very accurately asymptotic. It is somewhat daunting to see how com- 
manding the part 2tt plays as the controlling factor in the error for n. These circumstances 
have motivated our continuing inquiry into the analytic nature of the error function. 

An exact integral representation for the error, E(k), can be obtained. Returning to 
Eq. (128p for Rz^k), and following all of the procedure thereafter, we find 



8 



TABLE I: Values of E(k). 



k 


Using seven terms in Eq. (]35ap 


exact result 


5 


1.394929580622094159D-05 


1.458979572496761765D-05 


10 


1.208533257871241603D-10 


1.208426564955768191D-10 


15 


-1.717906040975859856D-15 


-1.717926356677219480D-15 


20 


-5.757361768118465909D-20 


-5.757368219957248525D-20 


25 


-5.054992959778042247D-25 


-5.054994756412840440D-25 


30 


6.170826902775628923D-30 


6.170826687049165966D-30 


35 


2.283000783226568157D-34 


2.283000802727741913D-34 


40 


2.2087769745461361 19D-39 


2.208776992026544785D-39 


45 


-2.093726806383270815D-44 


-2.093726801667894352D-44 


50 


-9.003884490090622337D-49 


-9.003884494714834077D-49 


55 


-9.605682948187048732D-54 


-9.605682954606388173D-54 


60 


6.756128797397263314D-59 


6.756128794821975275D-59 


65 


3.530031377831535807D-63 


3.530031378066361048D-63 


70 


4.137108141343023278D-68 


4.137108141778060533D-68 


75 


-2.030602678726272460D-73 


-2.030602678498624493D-73 


80 


-1.375355163972183061D-77 


-1.375355163990000725D-77 


85 


-1.764261661388176988D-82 


-1.764261661430246816D-82 


90 


5.380737425072894637D-88 


5.380737424803436715D-88 


95 


5.323448270450853269D-92 


5.323448270468033075D-92 


100 


7.454268189222908298D-97 


7.454268189274244908D-97 



Then, using the identity for the gamma function 



r(fc + 2) 

z k+2 



i 

oo 



t k+l e~ zt dL 



and interchanging orders of integration, gives 

Similarly, and following the exact same procedure, we have an integral representation for 
D2(k), and thus our final formula for the error, 

r °° cos(2ttx) x k e- ln ^ x §^dx 

E(k) - 8n i n (2)x r(x+i) T " ( 38 ) 

j x e v{x+\) ax 

To make the connection with this representation for E(k) and that given in Eq. (I35ap we 
observe that, for large values of k, the major contribution to the integrals comes from large 



9 



values of x. Therefore, Stirling's asymptotic expansion can be used to give 



T(x+ 1) 

r(x + |) 



I — 1 

+ 8x + I28x 2 



21 



1024x 3 32768x 4 



+ 



3=0 



(39) 



Inserting this expansion into Eq. ( 138|) . we immediately recover the series expansion for E(k) 
given in Eq. (I35ap where now Fj = 2 J A,-. To see clearly what this means, we turn to 
the valuable paper of Tricomi and Erdelyi [l9[ (see also Fields j20j) in which is given the 
asymptotic formula for the ratio of gamma functions, namely 



T(z + a) 



E 



r(* + /3) s~s j 



(_i)ir(^_ a + i ) +1 

T(/3 - a) J 



\arg(z + a)\ < n 



B (a 



a) x {a - p - j) ., 
0, 



(40) 
(41) 



where the B 



(a-/3+l), 



als j21|, defined by 



a) are the generalised Bernoulli polynomials, the Norlund polynomi- 



e* - 1 



3=0 



X 



\t\ < 2tt, 



(42) 



and where 



Bf\a-x) 



The coefficients, A,-, and hence the B f (x), can be calculated from the formula, Eq. (11) in 
ref. [19(| , or they can be computed numerically (as we have done). Using a = 1 and (3 
we find 



2 ' 



,(§) 



A, 



(2j - 3)\\ (-iy B^ (1) _ (2j-3)\\B 



2ij\ 2ij\ 
and comparing the Cj in Eq. (I25p with the B^ (x) in Eq. (1421) gives 



Therefore, Fj = 2^Aj is the same as —(2j — 3)!! Cj. The connection between the E(k) in the 
previous section and that herein is complete. 

Tricomi and Erderlyi opened their paper with the salient comment "Many problems in 
mathematical analysis require a knowledge of the asymptotic behavior of the quotient ffe^ 
for large values of There are two, in particular, that bear a kinship with our work. 

Watson 22| in a lovely cameo entitled A Note on Gamma Functions, studied tight bounds 
on the Wallis formula for it by employing the formula for the hypergeometric function, 
namely 



r(x + 1) 
r(x + i) 



1 2 



11 

x F(--,--;x;l) 



x 



(-1) (-5) 

V 2 / m V 2 / m 

J m(x + l) m -\ 



10 



with the comment, "the condition x + \ > amply secures the convergence of the series". 

In the classic problem of the one-dimensional random walk 23], the expectation value of 
the absolute distance after N (unit) steps is given by 



< dj > 



2 r(j + i; 
r(j) 



where J = ~ for N even, and 



N+l 



for N odd. With the asymptotic expansion for the ratio 



of these gamma functions, we recover the known result 



< d N >= \ j 



2N 



71 



T AN + 32N 2 



21 



128N 3 2048N 4 



where the top signs are taken for N even and the bottom ones for iV odd. The expansion for 
iV odd is precisely the one we have for the | expansion for the error, E(k). The alternating 
signs in < > is a direct consequence of the (— 1) J phase in Eq. (|42|) . We find the apparent, 
serendipitous, coincidence {Are There Coincidences In Mathematics? [24]]) quite intriguing. 



VI. EPILOGUE 



D. H. Lehmer was an eminent mathematician [25j. His many and varied works dealt 
especially with matters of numbers, the queen of mathematics, for which he had an abiding 
and prodigious affection and talent. Many of his elegant works bear his name in the litera- 
ture. Particularly relevant here are his Machin-like " On Arccotangent Relations for pi" 26] 
and " A Cotangent Analogue of Continued Fractions" 27|, in which he showed that every 
positive irrational number has a unique infinite continued cotangent representation. 

It would be remiss to close our essay without a salute to n, the most renown of all 
constants in mathematics. Of the literally plethora of fine possible choices, we have chosen 
to go with, 

" The value of tt has engaged the attention of many mathematicians and calculators from 
the time of Archimedes to the present day, and has been computed from so many different 
formulae, that a complete account of its calculation would almost amount to a history of 
mathematics. " 

[J. W. L. Glaisher, Messenger of Math., 25-30 (1872)] 

and 

" And he made a molten sea, ten cubits from the one brim to the other; it was round all 
about ... and a line of thirty cubits did compass it round about." 

[1 Kings 7:23] 

[On The Rabbinical Approximation of Pi 28] 
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TABLE II: The first 40 values of R 1 {k)/R 2 {k). 



1 


3 


0000000000000000000000000000000000000000000000000000000000000000, 


2 


3 


I42857l42857l42857l42857l42857l42857l42857l42857l42857l42857l429, 


3 


3 


1428571428571428571428571428571428571428571428571428571428571429, 


4 


3 


1415929203539823008849557522123893805309734513274336283185840708, 


5 


3 


1415780637940682708449916060436485730274202574146614437604924454, 


6 


3 


1415923909792426870850032900640066997667045522521983609499312078, 


7 


3 


1415928223305272485600354452813469206911829862649534780682321666, 


8 


3 


1415926602118958149513766625165465875475594263879799980191580879, 


9 


3 


1415926517076394721680972516575553164503064392560154252064551018, 


10 


3 


1415926534689505819670665641594107024726895639269814592980636856, 


11 


3 


1415926536097034333101766975356113602136332924145393370409221625, 


12 


3 


1415926535916966575323999453816550421791069920332492792677268321, 


13 


3 


1415926535895974772851022724935990443197089788766137673243522672, 


14 


3 


1415926535897658427283911681357586647531232166551257736863740840, 


15 


3 


1415926535897949563890000604989833119363870525182191538675275264, 


16 


3 


1415926535897936078697175611343357401374657772051447922473642141, 


17 


3 


1415926535897932264519119870605439217503823038833322854025505048, 


18 


3 


1415926535897932337393340424095495278361619490164656674686609330, 


19 


3 


-i A -i r AAAr arnAhTAnAfin ,( nr nnf?i r-^ r\ r\ c\ r" A ^ c\ n r\ r~t r~t c\ tr\ A r—r a t~\ r* c~i r-^ A c\ -t r—r-i rnA/inr r A/if? i 

1415926535897932384958871798265432629332247495874217158068559674, 


20 


3 


1415926535897932385202170654790753694551295971907030363946530365, 


21 


3 


1415926535897932384636042752252741477211700998636696041572103108, 


22 


3 


1415926535897932384619738974298389685253003898941919111583795322, 


23 


3 


1415926535897932384626159605817069391563876308322608250338128782, 


24 


3 


1415926535897932384626507820001679085203349744076397007956812907, 


25 


3 


1415926535897932384626438887789785254812134010734664381353470044, 


26 


3 


1415926535897932384626433063420855784182255379975982739476954955, 


27 


3 


1415926535897932384626433753704220177671845520815679347756243018, 


28 


3 


1415926535897932384626433840161803084753395271588726611787867400, 


29 


3 


1415926535897932384626433833921567846462444298200187678535709625, 


30 


3 


1415926535897932384626433832733320575101202334085977160126045251, 


31 


3 


1415926535897932384626433832780011313583089291827528827676225744, 


32 


3 


1415926535897932384626433832795412276163249733207758655099465346, 


33 


3 


1415926535897932384626433832795218564970456013498279015329014616, 


^4 


Q 
o 


1 41 l ^Q9fi c i' : l r i8Q7Q^9' : l84fi9fi4'^8' : i97Q^n9R7fi4Q1 9fi c lRfi71 fi4RD9 c > c i c i94fi c >9 c iRQXQ 


35 


3 


14159265358979323846264338327950265589708912660091443561059701 1 1 , 


36 


3 


1415926535897932384626433832795028788309975635209384389235382620, 


37 


3 


1415926535897932384626433832795028868138180463022120372825783062, 


38 


3 


1415926535897932384626433832795028843242936641656711090281944977, 


39 


3 


1415926535897932384626433832795028841687499915165439348592159580, 


40 


3 


1415926535897932384626433832795028841949606223830792761898241293, 
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TABLE III: Table M continued. The 41 st to 65 th values. 



41 3.1415926535897932384626433832795028841974582207305475602005895270 

42 3.1415926535897932384626433832795028841972028735133894452619084772 

43 3.1415926535897932384626433832795028841971667284657301213102341152 

44 3.1415926535897932384626433832795028841971689329076588332348800482 

45 3.1415926535897932384626433832795028841971694203123738376538881191 

46 3.1415926535897932384626433832795028841971694054855024473294830381 

47 3.1415926535897932384626433832795028841971693992698167077733580624 

48 3. 1415926535897932384626433832795028841971693992991 125354048682206 

49 3.1415926535897932384626433832795028841971693993745498910645297440 

50 3.1415926535897932384626433832795028841971693993760062094244160757 

51 3.1415926535897932384626433832795028841971693993751340966582560562 

52 3.1415926535897932384626433832795028841971693993750956731257005236 

53 3.1415926535897932384626433832795028841971693993751052376002021699 

54 3.1415926535897932384626433832795028841971693993751059289957208672 

55 3.1415926535897932384626433832795028841971693993751058305806275469 

56 3.1415926535897932384626433832795028841971693993751058199065942028 

57 3.1415926535897932384626433832795028841971693993751058208338153750 

58 3.1415926535897932384626433832795028841971693993751058209844036282 

59 3.1415926535897932384626433832795028841971693993751058209768684753 

60 3.1415926535897932384626433832795028841971693993751058209748770310 

61 3.1415926535897932384626433832795028841971693993751058209749198387 

62 3.1415926535897932384626433832795028841971693993751058209749448105 

63 3.1415926535897932384626433832795028841971693993751058209749448952 

64 3.1415926535897932384626433832795028841971693993751058209749445970 

65 3.1415926535897932384626433832795028841971693993751058209749445888 

7T 3.1415926535897932384626433832795028841971693993751058209749445923 



[1] D. H. Lehmer, AMM 92, 449 (1985). 

[2] L. Comtet, Advanced Combinatorics (Reidel, Dordrecht, 1974). 

[3] J. Borwein and R. Girgensohn, Aequat. Math. 70, 25 (2005). 

[4] R. J. Mathar (2009), arXiv:0905.0215vl [math.CA]. 

[5] T. Sherman (2000), http://math.arizona.edu/ rta/001/sherman.travis/series.pdf. 

[6] N. Batir, Appl. Math. Comp. 147, 645 (2004). 

[7] A. I. Davydychev and M. Y. Kalmykov, Nucl. Phys. B690, 3 (2004). 

[8] B. Sury, T. Wang, and F.-Z. Zhao, J. Integer Sequences 7, 2 (2004). 

[9] N. Batir, Proc. Indian. Acad. Sci. Math. Sci. 115, 371 (2005). 

[10] A. Sofo, J. Integer Sequences 9, 4 (2006). 

[11] R. Sprugnoli, Integers: Electronic J. of Combinatorial Number Theory p. 1 (2006). 

[12] J.-H. Yang and F.-Z. Zhao, J. Integer Sequences 9, 4 (2006). 



13 



TABLE IV: S k {2) for k = 1, . . . , 30 



1 3 + 7T 

2 11 + (7 7r)/2 

3 55 + (35 tt)/2 

4 355 + 113 7T 

5 2807 + (1787 vr)/2 

6 26259 + (16717 vr)/2 

7 283623 + 90280 vr 

8 3473315 + (2211181 vr)/2 

9 47552791 + (30273047 vr)/2 

10 719718067 + 229093376 vr 

11 11932268231 + (7596317885 vr)/2 

12 215053088835 + (136907048461 vr)/2 

13 4186305575415 + 1332542451241 vr 

14 87534887434835 + (55726440112987 vr)/2 

15 17/2 (230197719678574 + 73274209950655 vr) 

16 46561960552921315 + 14821132364094533 vr 

17 1175204650272267479 + (748158516941653967 vr)/2 

18 31357650670190565363 + (19962900431638852297 vr)/2 

19 881958890078887314567 + 280736233919792968780 vr 

20 26078499305918584929155 + (16602088291822017588121 vr)/2 

21 808742391638178302137783 + (514861397268710391722627 vr)/2 
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24 31443867356631172742458654755 + (20017787678934958873836057001 vr)/2 
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27 29/2 (120960271409361525575952065166694 + 38502850225074296050584562759015 vr) 

28 72107782245849606090651464405624515 + 22952619959641949809923983081211353 vr 

29 3068555154632012211023759300893608311 + (1953502884039199982724152729710551347 vr)/2 

30 135010171084427194623890031993168567507 + (85950144383076253408132013000868398677 vr)/2 



[13] J.-H. Yang and F.-Z. Zhao, J. Integer Sequences 10, 8 (2007). 
[14] A. Sofo, J. Integer Sequences 13, 2 (2010). 

[15] E. D. Krupnikov and K. S. Kolbig (1996), CERN/ Computing and Networks Document 
CN/96/8. 

[16] L. Euler, Institutiones calculi differentialis cum eius usu in analysi finitorum ac doctrina 
serierum, CAPUT VII (1775), reprinted from the collection of the University of Michigan 
Library. 

[17] A. P. Prudnikov et al., Integrals and Series, Vol.III (Nauka, Moscow, 1986), [see Eq. 7.3.6(1) 
and Table 7.3.1]. 

[18] N. D. Baruah, B. C. Berndt, and H. H. Chan, ANM 116, 567 (2009). 



14 



[19] F. G. Tricomi and A. E. Erdelyi, Pacific J. Math. 1, 133 (1951). 

[20] J. L. Fields, Proc. Edinburgh Math. Soc. 15, 43 (1966). 

[21] N. E. Norlund, Vorlesungen iiber Differenzenerechnung (Springer, Berlin, 1924). 

[22] G. N. Watson, Edinburgh Math. Notes 42, 7 (1959). 

[23] http://mathworld.wolfram.com/GammaFunction.html — for a nice coincise summary. 

[24] P. J. Davis, AMM 88, 311 (1981). 

[25] J. Brillhart, Acta Arithmetica, LX II.3 (1992). 

[26] D. H. Lehmer, AMM 45, 657 (1938). 

[27] D. H. Lehmer, Duke Math. J. 4, 323 (1938). 

[28] T. Boaz and D. Garber, Historia Mathematica 25, 75 (1998), 
http://u. cs.biu.ac.il/ tsaban/hebrew.html. 



15 



